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Abstract 

The phenomenology of a QCD-Pomeron model based on the ex- 
change of a pair of non-perturbative gluons, i.e. gluon fields with 
a finite correlation length in the vacuum, is studied in comparison 
with the phenomenology of QCD chiral symmetry breaking, based 
on non-perturbative solutions of Schwinger-Dyson equations for the 
quark propagator including these non-perturbative gluon effects. We 
show that these models are incompatible, and point out some possibles 
origins of this problem. 
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1 Introduction 



The interpretation of the Pomeron in the framework of QCD is not fully un- 
derstood. It is expected to be generated at least by two gluons exchange [||]. 
However, the exchange of two perturbative gluons cannot reproduce the ex- 
perimental results on diffractive scattering. In order to circumvent such 
difficulty, Landshoff and Nachtmann [0(LN) proposed a model where the 
Pomeron is described by the exchange of two non-perturbative gluons, whose 
properties are dictated by the expected structure of the QCD vacuum. These 
non-perturbative gluons should not propagate over long distances, i.e. there 
is a finite correlation length for the gluon field in the vacuum, which should be 
determined from first principles, and can be understood in terms of gluon con- 
densates [0. The LN model describes successfully the experimental results of 
diffractive phenomena using a quite simple ansatz for the non-perturbative 
gluon propagator [0, ^. Recently we improved the model using a QCD mo- 
tivated non-perturbative gluon propagator, where the correlation length is 
provided by a dynamically generated gluon mass 

The remarkable success of this Pomeron model might provide physical 
insight in the properties of the gluon propagator in the infrared region. It 
is, therefore, natural to speculate on the ability of such a model to describe 
other low-energy QCD phenomena, in particular, dynamical chiral symmetry 
breaking (DCSB). Within the spirit of the LN model, we study the com- 
patibility of both the Pomeron and DCSB phenomenologies. The quark 
self-energy is computed through the numerical calculation of the Schwinger- 
Dyson equation (SDE) for the quark propagator, employing the same non- 
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perturbative gluon propagators as the ones used in the Pomeron phenomenol- 
ogy. According to the traditional approach to the DCSB these self- 
energies are used to obtain the experimental values of the quantities that 
characterize the chiral breaking, such as the quark condensate and the pion 
decay constant. 

It is worth mentioning that the existence of a finite correlation length 
for the gluon field was already proposed ten years before the LN model in 
a completely different context 0, and it gained strong support by recent 
lattice simulation of the gluon propagator p], Therefore, a study of its 
phenomenological implications is quite compelling, although there are few 
places where this phenomenon can be effectively tested. 

Using only QCD motivated gluon propagators, i.e. propagators deter- 
mined from the solutions of Schwinger-Dyson equations or obtained from 
lattice simulation, which satisfy the LN conditions and are constrained by 
the diffractive scattering data, we show that they do not lead to satisfactory 
DCSB parameters. In Section II we discuss some of the non-perturbative 
propagators found in the literature, and show how they are constrained by 
the diffractive scattering data according to the LN model. In Section III we 
use the non-perturbative gluon propagators, with the parameters found in 
the previous section, to compute the relevant quantities for DCSB. In Section 
IV we discuss the incompatibility found in the previous section and point out 
the possibles origins of such a failure. 
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2 DifFractive scattering and gluon propaga- 
tors 



There are several experimental features of diffractive scattering that can be 
computed through the two-gluon exchange model of the Pomeron according 
to the LN prescription, and here we will choose two simple quantities which 
can be related to hadron-hadron total cross sections, and show how they 
constrain the gluon propagators. The first one will be the quark-Pomeron 
coupling and the second one, the pion-proton total cross section. 

The strength of the Pomeron coupling to quarks (/5o) at leading order is 
given by 



367r2 



(1) 



where g'^{K'^)/ATT is the running quark-gluon coupling, and the value of (3^ = 
4 GeV~^ is extracted from the proton-proton cross section. Notice that we 
will be working with Euclidean momenta {K"^ = — A;^). 
The amplitude of meson-nucleon scattering is |]TD 
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A = i—sfd^K[aD{K^)][aD{{2Q-Kf)] 
9 J 

x2[/M(g2)-/M((Q-K)2)] 



x3 



fNiQ')-fNQ'-T;Q-^ + ^K 
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(2) 



where s is the square of the center-of-mass energy, to the couplings (a) we 
associate the same momentum of its multiplying propagator, /m and /at are 
respectively the meson and nucleon form factors. The total cross section is 
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related to this amplitude by 

ImA{s,t = 0) 

or = , (3j 

s 

and Eq.(|D is computed with the form factors in the pole approximation 

^(^') = / X - (4) 



We will calculate the pion-proton total cross section. Actually, the two 
gluon exchange model cannot describe the full cross section growth with 
energy (oc sO-o^os-j^ which appears in a recent fit for the tt — p total cross 
section [IT; 



irp 



The model can only accommodate the value 13.63 in the above expression. 
We used the following mean squared radii for the proton and pion, respec- 
tively IH: < >= 0.67 fm^ and < >= 0.44 fm^. Inspection of Eq.(||) 
and Eq.(0) show that they are quite dependent on the gluon propagator ex- 
pression at A;^ — s> and, therefore, well suited to study its infrared behavior. 

We will now consider the propagators that satisfy some of the basic as- 
sumptions of the LN model: (a) the propagator has a finite correlation length 
and (b) it is finite at k"^ = 0. We will only deal with QCD motivated propa- 
gators, i. e. those obtained as solutions of the SDE for the gluon polarization 
tensor, as well as obtained in a lattice simulation. We start with the prop- 
agator determined by Cornwall [0, which has already been used by one of 
us to describe some diffractive scattering processes M , and has the following 
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expression 



A2 



(6) 



where im?{K'^) is a momentum-dependent dynamical mass 



In^ 



-12/11 



(7) 



In the above equation rrig is the effective gluon mass, and b = [UN — 
2nj)j^%T\'^ is the leading order coefficient of the (5 function of the renormal- 
ization group equation, with = 3 for QCD, and where n/ is the number of 
flavors taken to be 3. In Eq.(^), the quark-gluon coupling strength g is such 
that g ~ 1.5 — 2, which was determined from a fit of Eq(P) to the numeri- 
cal solution of the SDE for the gluon propagator. For the running coupling 
constant giiC^^ we assume a functional form which interpolates between a 
constant and the renormalization group asymptotic behavior 

g^{K'^) _ 127r/(33 - In^) 



(8) 



where r = k^/A^, and we chose k 



rrir, 



Such freeze-out of the coupling 



constant in the infrared region is consistent with the Cornwall and Papavas- 
siliou [|l^] study of the trilinear gluon vertex. A similar form for the coupling 
constant has also been used in many phenomenological applications (see, for 
instance, Ref. and references therein). 

It has also been pointed out that the dynamical gluon mass may have 
a faster decrease with the momentum [|l^, according to the operator prod- 
uct expansion (OPE) determination of the ultraviolet behavior of the gluon 
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polarization tensor |T6 

n (k'\ 34iV7r {a,G') 

where (asC^) is the gluon condensate. Therefore, to be consistent with the 
massive Cornwall propagator and with OPE, we will consider a gluon self- 
energy that interpolates between the constant infrared behavior of Eq.(|^) 
and the ultraviolet one of Eq.(|^) which will be given by 

m\K') = ^^^(/.^ - K') + ^^e{K' - (10) 

where the scale /i^ will be limited by Eq.(|l]) and Eq.(0). 

Another infrared finite propagator has been found by Stingl and collab- 
orators as a solution of the SDE for the gluon polarization tensor fl^. Its 
form agrees with that derived by Zwanziger based on considerations related 



to the Gribov horizon |18|, and is given by 



DJK^) = (11) 



K2 

,4 / 1^2 



where /i^ is a mass scale not determined in Ref. ^T^. The term fig/K is 
exactly what is expected by OPE analysis whenever a mass scale for the 
gluon is introduced, and the full solution found in Ref. did contain a 
mass term, although its consequences were not pursued. It is also interesting 
to notice that the Bernard et al. lattice result for the gluon propagator 
can be fitted by Eq.(||) as well as Eq.(|lT]). 

Finally, Marenzoni et al. P] also performed a lattice study of the gluon 
propagator in the Landau gauge, obtaining for its infrared behavior the fol- 
lowing fit 

Dn,{K^) = 77^, (12) 
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where mi, Z and t] are constants determined by the numerical simulation. 
mi is of 0{A ^ IQOMeV), Z ^ 0.4 and rj ^ 0.5 what is slightly different 
from the previous propagators. The results of Bernard et al. also show the 
behavior {K'^)^, but with a smaller value for 77. In this case all the propagator 
parameters are determined, and we simply have to see if they are consistent 
with diffractive scattering and DCSB. All the above propagators where ob- 
tained in Landau gauge, except Cornwall's one, whose massive solution was 
shown to be gauge invariant. 

Performing the calculation of the integrals in Eqs. (|l]) and (H) for each 
one of the above propagators, we obtained the curves for /3q displayed in 
Fig.(l), and the pion-proton total cross section shown in Fig. (2) as functions 
of the gluon mass. From these figures and the experimental values of /3o, and 
the pion-proton total cross section, it was possible to establish an effective 
gluon mass for each propagator. Our results are shown in Table ^ 

It is clear from this table that almost all the gluon mass scales that fit 
the experimental value of /3o give reasonable (within 10%) values for the to- 
tal pion-proton cross section, except in the case of the lattice propagator. 



Eq.(12), where the comparison with the experimental values yields effective 
gluon masses differing by 35%. It is evident that not only the mass scale 
plays a role in this phenomenology, but also the functional form of the prop- 
agator introduces differences in the calculation, and this is why this problem 
becomes more interesting, because we can speculate about the infrared be- 
havior of the gluon propagator. As we shall see in the next section, we will 
fail to obtain DCSB with the above propagators for the gluon mass scales 
presented in Table |I[ 



3 DCSB with LN type propagators 



We discuss dynamical chiral symmetry breaking following the traditional 
approach |p. This consists in solving the Schwinger-Dyson equation (SDE) 
for the quark propagator and look for a mass term dynamically generated in 
this propagator. 

The Schwinger-Dyson equation for the quark propagator is 

S-\p) =^-4/ ^,l,S{q)T^{p,q)g'D^%p-q), (13) 
where we write the gluon propagator in the form 

^ '-^^^ . ^) . (14) 

The propagator has been written in the Landau gauge, which will be used 
throughout our work. In the above equations T,^{p, q) is the vertex function, 
and a(— g^/A^) is the QCD running coupling constant given by Eq.(|]). 

To proceed further we also need to introduce an ansatz for the quark-gluon 
vertex T^{p,q), which must satisfy a Slavnov- Taylor identity that, when we 
neglect ghosts, reads 

{v-q),T'{v,q) = S-\p)-S-\q). (15) 

This identity constrains the longitudinal part of the vertex, and if we write 
S^^{p) in terms of scalar functions 

S-\p)=A{p)^-B{p), (16) 

we find the solution [|T^ 

r^(P,g) = ^^^,[W)-A{q')]4-[B{p')-B{q^)]) 

+A{p'^)'~^^ + transverse part, (17) 
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which is a much better approximation than the use of the bare vertex . 
Assuming that the transverse vertex part vanishes in the Landau gauge we 
obtain 

D^^'iv - g)r.(g,p) = D^\v - q)A{q^)lu, (18) 
and arrive at the approximate Schwinger- Dyson equation 

Wp^) - 11 ^ - Bip^) ^ ,\ I ^J^D^ip - (19) 

Going to Euchdean space, we will be working with the following nonlinear 
coupled integral equations for the quark wave-function renormalization and 
self-energy 

[A(p2)_i]p2 ^ IQt: f d^Q a{{P - Qf / 



3 J (27r)4 $[(p-g)2] 

P.{P-Q)Q.[P-QY 



X P.Q + 2 



{P-Q? 



(20) 



^ > J (27r)4 $[(P-g)2] A2(Q2)Q2 + p2(g2)' ^ ) 

where = —q^ and P^ = —p"^, and we introduced a function — Q)^] 

which, in the case of the perturbative propagator, is simply $[(P — Q)"^] = 
{P — QY, for a massive bare gluon it will have the form — Q)"^] = (P — 
+ and will be more complex expression according to the propagators 
we discussed in the previous section. 

The numerical code we used to solve the above equations is the same of 
Ref. ||21[|, and for each one of the propagators in Section II we substitute: 
(a) ^K') = D-\K% (b) ^K') = D-^iK'), (c) ^K') = D;\K') and 
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(d) = D^i^K"^). After obtaining for each case the function A(p^) and 

B{j>^), we can compute the quark condensate, which is expressed as 

where = A{p'^) and M{p'^) = B{p'^)/A{p'^), and we will also calculate 

the pion decay constant, which in terms of the functions A and B of Eq. (|20|) 
and Eq.(pTD is given by 

f2 - 10, / ( p^d\nA\ 

p'^f^dA JB\(^ 2^^1n^\n 

For all the propagators discussed here, with the respective gluon mass 
scales given in Table |I| , we have not obtained dynamical quark mass gen- 
eration! Our results for f^^ and {qq) are identically zero (at least ten orders 
of magnitude below the scale A). The gluon mass scales of Table |l| are too 
large, in the sense that they cause a too strong screening of the force neces- 
sary to generate the symmetry breaking. To illustrate this incompatibility 
we present in Fig. (3) the curve of the dynamically generated quark mass 
(M(p^ = 0)) against the gluon mass, in the case of the Cornwall propagator 
(Eq.(ID), which gives the largest signal of mass generation that we obtained 
among all the propagators described above. Notice that in Fig. (3) to obtain 
dynamical quark masses of (9(300 MeV), we would need a gluon mass almost 
half the value necessary to satisfy the constraints from diffractive scattering. 

The inconsistency between diffractive scattering and the chiral symmetry 
breaking phenomenologies is not only a matter of adjustment of the infrared 
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gluon mass scale, this scale really plays a different role in both cases. We 
need large gluon masses {rrig » A) to soften the t-channel singularity in 
diffractive scattering, but smaller gluon masses if we do not want to erase 
the dynamical quark mass generation. It is obvious that we have an incon- 
sistency, although it is far from clear which is the solution. Nevertheless, 
this comparison of phenomenologies gives a powerful tool to constrain the 
detailed behavior of the infrared gluon propagator. 

4 Conclusions 

The LN Pomeron model is able to explain a large amount of experimental 
data on diffractive scattering 0, |^, making use of an infrared finite gluon 
propagator. As we have seen in Section 2, we need only a unique gluon mass 
scale to fit the pomeron-quark coupling and the pion-proton total cross sec- 
tion. On the other hand we have an extensive and successful phenomenology 
of DCSB, which was performed mostly with a "perturbative" gluon 
propagator. In Section 3 we followed the standard procedures to compute 
the dynamical quark mass, quark condensate and pion decay constant, with 
the same non-perturbative propagators prescribed by the LN model. Both 
phenomenologies are inconsistent. We need large gluon masses for diffractive 
scattering, and small ones to generate appropriate chiral symmetry breaking 
parameters! 

There is a clear difference between our calculation of DCSB with the 
previous ones. The use of "massive" propagators clearly screens the force 
necessary to generate the chiral symmetry breaking, and it is probably in 
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this direction that we may look for troubles, another possibility would imply 
that the LN model is incorrect, or that the standard phenomelogy of DCSB 
has nothing to do with a finite correlation length for gluons. Let us dis- 
cuss about each one of these possibilities. The first common point to both 
phenomenologies is the existence of a mass scale for the gluon propagator, 
however, there is another claimed form for the gluon propagator in the in- 
frared region, namely, (see |22| and references therein). In view of the 
theoretical arguments of Ref. [0, |16|, [l^, |l8l about gluon mass scales, and 



the numerical simulations on the lattice ^, we were compelled to assume 
that the gluon field definitively has a finite correlation length in the vacuum, 
discarding the solution. 

The gluon propagator being infrared finite, there is no comparison with 
the experimental data which would force us to abandon the LN model, and 
it remains to see if the DCSB phenomenology with such propagator is still 
consistent. Another difference to standard DCSB calculations is the use of 
the QCD running coupling constant given by Eq.(^, which, as we said pre- 
viously, is consistent with a study of the trilinear gluon vertex, and naturally 
embodies the freezing of the coupling at the relevant gluon mass scale [1^ 



Both phenomenologies, diffractive scattering and DCSB, depend strongly on 
the behavior of this coupling in the infrared, however, its effect works basi- 
cally in the same direction, i.e. increasing the coupling in the infrared we 
increase the diffractive scattering parameters as well as the dynamical quark 
mass. Therefore, there will not be much room to changes in this direction. 
Finally, if we assume that the LN Pomeron model is correct, the gluon field 
has a finite correlation length in the vacuum, and the coupling constant has 
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the freezing-out referred to above, we are forced to say that the dynamical 
quark mass is not generated by single-gluon forces, and the usual calculation 
of dynamical quark masses should be modified when the gluon mass is taken 
into account. Actually, this view is not new, studying the coupled fermion 
gap and vertex equations for DCSB in QCD it was also found an incompati- 
bility pSf , and it was argued that the standard technique of using one-gluon 
exchange in the SDE is not suitable when we consider the effect of massive 
gluons, and our work, in a different way, corroborates their result. 

We finalize stressing that we are still far from knowing all the sub- 
tleties of QCD, and we have formalized the incompatibility between a QCD 
Pomeron model and DCSB. We believe that the confront of these different 
phenomenologies may provide strong constraints on the infrared behavior of 
the gluon propagator. 
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Figure Caption 

Table 1 Gluon masses (in MeV) obtained fitting f3o and aj^p for tlie several 
propagators discussed above. We denote by Dcm the gluon propagator given 
by Eq.(|^) with the dynamical gluon mass given by Eq.([TD|). 

Fig. 1 Curves for (3q as a function of the gluon mass. We constrain the 
gluon mass by fixing (3^ = 4.0GeV~^. Each gluon propagator is represented 
by: dashed line (-De), dash-dotted line (-Dm), dotted line (-D^), solid line 

(-Dmc) • 

Fig. 2 Curves for pion-proton total cross section as a function of the 
gluon mass. The constraint on the gluon mass comes from fixing the value 
(j^p = 13.63 mb. The gluon propagators are represented by: dashed hne 
(Dc), dash-dotted line (Dm), dotted line (Dmc), solid line (-D^). 

Fig. 3 Dynamically generated quark mass for the Cornwall's gluon propaga- 
tor function of the gluon mass. 
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Figure 1: 
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Figure 2: 
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Figure 3: 
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